Design of Wideband Sampled-Data Filters 

By R. M. GOLDEN and J. F. KAISER 

(Manuscript received March C, 196-1) 

A design procedure is presented for readily obtaining sampl ed -data ji Iter 
representations of continuous filters. The procedure utilizes the bilinear z 
transformation and preserves the essential amplitude characteristics of the 
continuous filter over the frequency range between zero and one-half the 
sampling frequency. It is shown that the procedure can yield meaningful 
sampled-data filter designs for many of those filters where the standard z 
transform cannot be used directly. 

I. INTRODUCTION 

Sampled-data filter representations for continuous filters can be ob- 
tained using several different design procedures. 1 A particular design 
method utilizing the bilinear transformation is developed herein. The 
method is especially useful in designing wideband* sampled-data filters 
which exhibit relatively flat frequency-magnitude characteristics in suc- 
cessive pass and stop bands, filters of this type are widely used in net- 
work simulation and data processing problems.'-' The design method pos- 
sesses two chief advantages over the standard z transform. 3 The first is 
that the transformation used is purely algebraic in form. This means it 
can be applied easily to a continuous filter having a rational transfer char- 
acteristic expressed in either polynomial or factored form. The second 
advantage is the elimination of aliasing 4 errors inherent in the standard 
z transform. Thus, the sampled-data filter obtained by this design 
method exhibits the same frequency response characteristics as the con- 
tinuous filter except for a nonlinear warping of the frequency scale. 
Compensation for this warping can be made by a suitable frequency 
scale modification. Some of the more common filter networks to which 
the design method can be applied effectively are the Butterworth, Bessel, 
Chcbyshcv, and elliptic filter structures. 

The essential properties of the bilinear transformation are presented 

* A sampled-data filter design will he termed "wideband" if the frequency range 
of useful approximation approaches half the sampling frequency. 
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in the next section. For comparison purposes, properties of the standard 
z transform are also given. This is followed by a detailed description 
of a filter design procedure using frequency transformations. Examples 
illustrating the design procedure arc then presented. A short discussion 
concerning computer simulation of the obtained sampled-data filters is 
also included. 

II. THE STANDARD AND BILINEAR 2 TRANSFORMATIONS 

In this section it is assumed that a satisfactory rational expression is 
known for the transfer function of a continuous filter for which a sampled- 
data approximation is sought. What is then necessary is a transformation 
which will convert this transfer characteristic into a sampled-data trans- 
fer function rational in z~ l , the unit delay operator. 

Two transformations applicable to this problem are the standard z 
transform and the bilinear z transformation. 

The standard z transform applied to H(s), the transfer function of the 
filter, is 3 

H*(s) = E H(s+jmu> a ) (1) 

m=— oo 

or equivalently in terms of the impulse response, h(t), of the filter 



3C*(z) = TY,KlT)z- 1 (2) 

i=0 



where 



s = a + ju 
H(s) = Laplace transform of h(t) 

cos = 2ir/T = radian sampling frequency 
H*(s) = Laplace transform of the sampled filter impulse re- 
sponse 
z~ l = exp (-sT) = the unit delay operator 
3C*(z) = H*(s) | „ = (in*)/r = 2 transform of h(t). 

The behavior of H(s) for s greater than some critical frequency jo) e is 
assumed to be of the form 

H(s) 1.1. „> ,«„ = K/s n , n > (3) 

where K is a determined constant. 

Equation (1) or (2) is the transfer function of a sampled-data filtei 
which approximates the continuous filter. In the time domain, the im- 
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pulse response of the sampled-data filter is the sampled impulse response 
of the continuous filter. This can be shown by taking the inverse trans- 
form of (2). Equation (1) shows that in the baseband 

(-W./2 ^ u ^ u«/2) 

the frequency response characteristics of the sampled-data filter, H*(s), 
differ from those of the continuous filter, H(s). The difference is the 
amount added or "aliased" 4 in through terms of the form 

H(s + jmut), m 9^ 0. 

If H(s) is bandlimited to the baseband, i.e., | H(s) | = for u> > a>„/2, 
then there is no aliasing error and the sampled-data filter frequency 
response is identical to that of the continuous filter. Unfortunately, 
when H(s) is a rational function of s, it is not bandlimited and therefore 
II (s) ?* H*(s) in the baseband. 

The magnitude of the errors resulting from aliasing is directly related 
to the high-frequency asymptotic behavior of H(s) as defined in (3). If 
n is large and u> c « w,/2, then the aliasing errors will be small and the 
standard z transform generally will vield a satisfactory sampled-data 
filter design. However, in wideband designs, w c is usually an appreciable 
fraction of u>„/2. Furthermore, many continuous filter designs result in 
transfer functions in which n is 110 greater than 1. These two conditions, 
namely co,. ~ w s /2 and n = 1 , can create large aliasing errors in the 
frequency response characteristics, thus yielding an unusable result. 

Fortunately, even when m c « a>„/2 and n = 1, a design method em- 
ploying the bilinear z transformation may provide satisfactory wide- 
band designs. This z form is defined from the mapping transformation, 

* = (2/T) tanh (siT/2) (4) 

where 

Si =(Ti+ jwi . 

The right-hand side of (4) is periodic in u>i with period 2-ir/T. Consider- 
ing only the principal values of wj , —w/T < «i < r/T, it is seen that 
the transformation given by (4) maps the entire complex s plane into 
the strip in the Si plane bounded by the lines on = —ir/T and wi = 
-\-ir/T. For this reason the bilinear transformation can be looked upon 
as a baiidlimiting transformation. Therefore, when this transformation 
is applied to a transfer function H(s), the entire s-plane frequency 
characteristics of H(s) are uniquely carried over into the frequency 
characteristics of H(si). 

* This transformation will be referred to as the bilinear 2 form or 2 form. 
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With the substitution 



-1 —liT 

z = e , 



(4) can be written immediately as 



s = 



2 (1 - z~ l ) 



(5) 



T (1 + z- 1 ) * 
Thus, 

0C(z) m H(s) I 2 ( i_ 2 -i) (6) 



r 



(1+z-i) 



where K(z) denotes a sampled-data transfer function obtained by the 
use of the bilinear z form.* 

The transfer function 0C(z) obtained by means of the bilinear z form 
and the function W*(z) obtained by means of the standard z transform 
are both rational in z~ l and of the same denominator order as the con- 
tinuous filter. These two functions, 3C(z) and 3C*(z), become essentially 
equal to each other as the sampling frequency becomes large compared to 
the moduli of each of the poles of the continuous filter function, H(s). 
When the sampling frequency is not large, representation of some filters 
by the standard z transform can be quite unsatisfactory because of 
aliasing errors. However, the bilinear z form, with its absence of aliasing 
errors, may give a satisfactory representation for these filters. A particu- 
lar set of filters to which the bilinear z form always can be applied suc- 
cessfully are those which exhibit relatively flat frequency-magnitude 
characteristics in successive pass and stop bands. This follows directly 
from (6). 

Thus, sampled-data filters designed by using the bilinear z form pre- 
serve the essential amplitude characteristics of the continuous filter. In 
the baseband ( -w s /2 ^ w ^ w,/2), the frequency characteristics of the 
sampled-data filter are identical to those of the continuous filter except 
for a nonlinear warping of the frequency scale. 

This warping is found from (4) upon substituting ju for s and is 

co = (2/T) tan (^T/2). (7) 

For small values of wi , the relation is essentially linear, producing 

* It should be noted that the bilinear transform is used here in a distinctly dif- 
ferent way than it is commonly used in sampled-data control system design. In 
the control system literature it is used to transform the sampled-data function 
JC*(2) from the discrete domain back to the continuous domain for conventional 
stability and frequency response analysis. See for example J. T. Tou, Digital and 
Sampled-Data Control Systems, McGraw-Hill, New York, 1959, pp. 244-247 and 
pp. 466-470. 
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Fig. 1 — The frequency scale warping of the bilinear z transformation. 

negligible warping at the lower end of the frequency scale. Fig. 1 shows 
the nature of this warping. Compensation for the effect of warping can 
be made by prewarping the band-edge frequencies of the continuous 
filter in such a way that application of the z-form transformation will 
shift the band-edge frequencies back to the desired values. The incor- 
poration of this prewarping compensation into a sampled-data filter 
design procedure is discussed in the next section. 



III. A SAMPLED-DATA FILTER DESIGN' METHOD 

A sampled-data filter design may be obtained by applying the 2-form 
transformation of (6) to the rational transfer characteristic for a continu- 
ous filter. However, in order to compensate for the frequency warping 
imposed by the z form, the frequency characteristics of the continuous 
filter first must be altered or prewarped. Hence the transfer characteristic 
for the continuous filter must be redesigned such that the band-edge 
(cutoff) and maximum loss frequencies are computed according to, 



co c = (2/7 1 ) tan (u d T/2) 



(8) 
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where: 

a) c = computed cutoff or loss frequency 
wd = desired cutoff or loss frequency. 

The redesign of the continuous transfer characteristic cannot be ac- 
complished simply by applying (8) to each pole and zero of the original 
transfer characteristic. On the contrary a completely new transfer char- 
acteristic must be obtained for the continuous filter. It is then possible 
to obtain the desired sampled-data filter by applying the 2-form of (6) 
directly to the redesigned transfer characteristic of the continuous filter. 
The sampled-data filter so obtained will then have the desired magnitude- 
frequency characteristics. 

Compensation for frequency warping becomes especially simple to 
apply if the original continuous filter design was obtained by applying 
a frequency-band transformation to a suitable low-pass design such as 
Butterworth, Chebyshev, etc. Thus the extensive literature available on 
tabulated low-pass filter designs can be used to great advantage to sim- 
plify the filter design problem. The well-known frequency transforma- 
tions which convert a normalized low-pass filter to a low-pass, a band- 
pass, a bandstop, or a high-pass design are 

s n = s/oiu low-pass to low-pass (9) 

s„ = , """! low-pass to bandpass (10) 

s(to u — on) 

s n = s } Uu - Ul \ low-pass to bandstop (11) 

(S- + WuWi) 

s n = 03,,/s low-pass to high-pass (12) 

where : 

s n = the complex variable of the normalized low-pass filter transfer 

function 

8 = the complex variable of the desired filter transfer function 
co„ = the upper radian cutoff frequency 
u>i = the lower radian cutoff frequency. 

When continuous filters are designed with the aid of these transforma- 
tions, prewarping is accomplished by properly choosing the cutoff fre- 
quencies used in the frequency transformations. The choice of these cut- 
off frequencies is determined from the desired cutoff frequencies by 
means of (8). Using these values, the new prewarped transfer f miction is 
determined by applying (9), (10), (11) or (12) to the original low-pass 
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function. Transformation is made to a sampled-data filter by applying 
(5) to the prewarped continuous function. This sampled-data filter will 
now have the correct cutoff frequencies. The transfer function thus ob- 
tained can be used directly in a digital computer simulation. 

IV. SIMULATION OF SAMPLED-DATA FILTERS 

Application of either the standard z transform or the bilinear z form 
to a rational transfer function yields a transfer function rational in z~ [ 
for the sampled-data filter. The programming or simulation of this 
sampled-data filter on a digital computer can be accomplished by either 
the direct, the cascade or the parallel form. These forms, as commonly 
defined, are shown in Fig. 2. In this figure 9(2) and 5(z) represent finite 
polynomials in z~ l for feed-forward and feedback transmissions re- 



ft) DIRECT FORM 



3(z) 



J 



9G0 



3C(z) = 



I + 2FU) 
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(C) PARALLEL FORM 
Fig. 2 — Some possible simulation forms for sampled-data filters. 



1540 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1964 

spectively; where subscripted, the order of the polynomial is at most 
second. 

The choice of which of the three forms to use for simulation of the 
.sampled-data filter depends on the complexity of the filter function 
H(s), on the form of H(s), and on the particular z transformation used. 
Generally, simulation by the direct form requires considerably greater 
accuracy in the determination of the filter parameters than either of the 
other two forms. This is especially true when the order of H(s) is large 
and when H(s) has poles with real parts that are a very small fraction 
of the sampling frequency. For this reason either the cascade or parallel 
form may be preferred. 

The choice between using the cascade or the parallel form depends 
largely on which 2-transform method is used to obtain the sampled filter 
and how that particular method is applied. Realization in the cascade 
form requires calculation of the numerator polynomial, Q(z), or its fac- 
tors. This computation consists of a simple algebraic substitution when 
the bilinear z form is applied to a filter function H (s) expressed in the 
form, 

G(s) 

11 - a k ) 
fe-i 

Determination of Q(z) in polynomial or product form respectively allows 
either of the following cascade realizations to be made: 



or 



_ £ / W + ouT 1 + a^ \ (15) 

£=1 \ 1 + hkZ - 1 + buz 2 / 

If the numerator G(s) is in polynomial form, considerable care must be 
taken in the calculation of the coefficients of the polynomial £(z), as this 
computation involves differencing nearly equal numbers. 

For realization in the parallel form the partial fraction expansion of 
H(s) must be known. Since in the standard z transform method obtain- 
ing the partial fraction expansion is a necessary step, simulation of 
filters designed by this method is most directly accomplished in the 
parallel form. Here the continuous filter transfer characteristic is repre- 
sented by 
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H(s) = £ 



£1 (hs 2 + Qu-s + Qo* 
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Transforming this expression by use of the standard z transformation 
yields 



K(z) = Z 



4,*2 + A, 



(17) 



£3 B 2t r- S + f* l/t2 -» + 1 

whereas transforming by use of the bilinear z form yields the similar 
expression 



X(z) = (1 + z" 1 ) E 



Au*~ x + Ao,- 



£{ B 2 ,z- 2 + Buz-' + 1 • 



(18) 



Each rational function in either of the above summations can be syn- 
thesized by the recursive structure shown in functional block diagram 
form in Fig. 3(a). This recursive structure uses only two delays, four 
multiplications, and live additions. The complete realization of (18) is 
shown in Fig. 3(b). 




3Cj(z) = 



A,iZ-« + A 



°j 



B 



, jZ -2 + B,jZ-<+1 



(a) 

TERM IN PARTIAL 
FRACTION EXPANSION 



3C(z) = + z-') Z JCj(z) 



(b) 

COMPLETE TRANSFER FUNCTION 

EXPRESSED AS A PARTIAL 

FRACTION EXPANSION 



Fig. 3 — Simulation in parallel form of ;i sampled-data filter obtained by the 
bilinear z transformation. 
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The programming of these sampled-data filters for computer simula- 
tion can be greatly simplified if a compiler such as the Block Diagram 
(BLODI) compiler 5 developed at Bell Telephone Laboratories is used. 
The compiler permits specification of a sampled-data system in func- 
tional block diagram form. 

In the following section an example is presented for a filter designed, 
synthesized, and simulated by the foregoing method. 

V. DESIGN, SYNTHESIS, AND SIMULATION OF A WIDEBAND BANDSTOP 
SAMPLED-DATA FILTER 

As an example of the application of the bilinear z-form to the simula- 
tion of a practical filter, consider the design of a particular bandstop 
filter. The filter is to exhibit at least 7") db loss in a rejection band which 
extends between 2596 cps and 2836 cps. Below 2588 cps and above 2844 
cps, the loss is to be between and +0.5 db. The sampling rate of the 
discrete filter is to be 10 kc. The complexity or order of the filter is to be 
held to a minimum. Fig. 4 shows a sketch of the amplitude response char- 
acteristics desired of the filter. 

Minimum filter complexity and sharp transition between pass and stop 
bands suggest the use of an elliptic filter 6 (equiripple) as the basic low- 
pass type. However, before a suitable low-pass structure can be deter- 
mined, the above specified critical frequencies must be prewarped by 




2.588 | | 2.844 

2.596 2836 

FREQUENCY IN KILOCYCLES PER SECOND 



Fig. 4 — Desired amplitude response characteristic of a bandstop filter. 
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means of (8) . The warped values are : 

f,„ = lower cutoff frequency in passband (2588 cps) = 3364.15 cps 

fir — lower cutoff fre([iiency in rejection band (2596 cps) 

= 3381.13 cps 

f ur = upper cutoff frequency in rejection band (2836 cps) 

= 3937.54 cps 

f up — upper cutoff frequency in passband (2844 cps) = 3957.84 cps. 

The warped values at the lower band edge require a low-pass filter 
with a transition ratio of 0.93792, while the values at the upper band 
edge require a transition ratio 6 of 0.93658. Therefore, to meet the original 
specifications, the larger of the two transition ratios must be chosen. 
Hence specifications required for the basic low-pass elliptic filter are: 

in-band ripple = 0.5 db 

out-of-band minimum attenuation = 75.0 db 

transition ratio = 0.93792. 

Application of elliptic filter design procedure with these specifications 
yields a basic low-pass structure of eleventh order. The poles and zeros 
for the transfer function of this low-pass filter are listed in Table I. The 
low-pass filter has been normalized to have a cutoff frequency of one 
radian per second and amplitude gain of unity at zero frequency. 

Table I — Poles and Zeros of Normalized Eleventh-Order 
Elliptic Low-Pass Filter 



In-band ripple = 0.500 db 
Minimum attenuation = 76.501 db 
Transition ratio = 0.937917 

Gain factor = 0.001 1060 

Poles 

-0.0069130 ±j 1.0010752 
-0.0257C.Ui ±j 0.9750431 
-0.0615122 ±j 0.9063786 
-0.1269215 ±j 0.7504391 
-0.2142976 ±j 0.4483675 
-0.2611853 

Zeros 

±j 1.0695414 
dbj 1.1009005 

±i 1.1946271 
±j 1.4652816 
±j 2.5031313 
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The desired bandstop filter is obtained next by applying the low-pass- 
to-bandstop transformation, given in (11), to the normalized elliptic 
low-pass filter. The cutoff frequencies used are the warped values ob- 
tained above for fi p (3364.15 cps) and f up (3957.84 cps). The resulting 
bandstop filter is then transformed by the bilinear z form to yield the 
required sampled-data filter. Table II lists the coefficients of the result- 
ing sampled-data filter needed for parallel realization of the form shown 
in Fig. 3. Fig. 5 shows the frequency response characteristics of this 
sampled-data filter. It is seen that the original filter specifications are 
met by the sampled-data filter. For comparison purposes, the standard 
z transform was applied directly to the twenty-second-order continuous 
filter. The frequency response characteristic of this filter is shown in Fig. 
6. It is seen that the standard z transform has yielded an unusable result. 

VI. SUMMARY 

The need for sampled-data filters in wideband simulations of many 
processing systems has led to a synthesis method which overcomes the 
shortcomings of the standard z transform. The method presented con- 
sists of directly transforming a suitable continuous transfer function to a 
sampled-data filter by means of the bilinear z form. For wideband filters 
the method is particularly suited to those filters that exhibit relatively 
constant magnitude-frequency characteristics in successive pass and 
stop bands. Conventional design techniques of continuous filters are used 



Table II — Partial-Fraction Expansion Coefficients for 

Parallel Realization of Sampled-Data 

Bandstop Filter 





Numerator Coefficients 


Denominator Coefficients 




Ai 


At 


Bi 


Bi 


1 


0.0001628 


0.0008827 


0.9987854 


0.1106416 


V 


-0.0009283 


-0.0001764 


0.9989898 


0.4285348 


3 


-0.0024098 


-0.0027894 


0.9956089 


0.1063723 


4 


0.0031774 


0.0026966 


0.9957459 


0.4317548 


ft 


0.0102446 


0.0026026 


0.9879911 


0.0940731 


fi 


-0.0037799 


-0.0112135 


0.9883051 


0.4414974 


7 


-0.0277640 


0.0127415 


0.9651789 


0.0616261 


8 


-0.0108027 


0.0289421 


0.9661438 


0.4663508 


q 


0.0272223 


-0.1163873 


0.8694592 


-0.0204564 


10 


0.1206914 


. -0.0054765 


0.8742300 


0.5186036 


n 


0.2973946 


-0.2973227 


0.5283651 


0.2074591 
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Fig. 5 — Frequency response characteristics of the sampled -data bandstop 
filter designed by the bilinear z transformation. 

directly in the synthesis procedure of the sampled-data filters. (Thus the 
synthesized filters have frequency characteristics comparable to those of 
continuous filters.) An example has been presented of a filter function 
synthesized by this procedure and easily programmed for a simulation. 
Results obtained from this example demonstrate the usefulness and ac- 
curacy of the bilinear g-form method. 
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Fig. 6 — Frequency response characteristics of the sampled-data bandstop 
filter designed by the standard z transformation. 
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